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The  practical  consequences  of  the  commonly  occurring  thermospheric  condition  of 
turning  points  for  propagating  acoustic  and  internal  gravity  wave  modes  have  been 
investigated.  This  has  required  some  rather  novel  theoretical  developments  on  the 
W.K.B.  approximation,  turning  points,  and  the  measurement  of  phase  velocities.  In 
the  process  of  carrying  out  this  investigation  some  sample  calculations  were  made 
illustrating  the  relationships  between  phase  velocity  measurements  and  the  existing 
values  of  the  Vaisala  frequency  and  acoustic  cut-offs.  A  study  was  also  made  of  the 
reflection  of  internal  gravity  waves  of  long  period  at  height  heights  for  which  the 
wavelength  becomes  of  the  order  of  the  mean  free  path  of  the  molecules.  The  net 
effect  is  one  of  damping  on  the  long-range  propagation  of  trapped  internal  waves. 
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i he  program  of  research  offered  in  the  work  proposal  for  F.Y.  1971 
(December  1,  1970  -  November  30,  1971)  has  been  fulfilled  as  follows: 

1)  We  have  investigated  the  practical  consequences  of  the  commonly  occur¬ 
ring  thermospheric  condition  of  turning  points  for  propagating  acoustic 
and  internal  gravity  wave  modes.  This  has  required  some  rather  novel  theo¬ 
retical  developments,  which  we  have  entitled  :  "The  W.K.B.  Approximation, 

Turning  Points  and  the  Measurement  of  Phase  Velocities"  (see  below) . 

♦ 

2)  In  the  process  of  carrying  out  1),  we  have  made  some  sample  calculations 
illustrating  the  relationships  between  phase  velocity  measurements  and  the 
existing  values  of  the  VaisSia  frequency  and  acoustic  cut-offs  (figures  1, 

2.  3). 

3)  We  have  also  included  a  very  brief  and  probably  oversimplified  study  of 
the  reflection  of  internal  gravity  waves  of  long  period  at  great  heights 
for  which  the  wavelength  becomes  of  the  order  of  the  mean  free  path  of  the 
molecules.  Since  the  net  effect  is  one  of  damping  on  the  long-range  propa¬ 
gation  of  trapped  internal  waves,  we  have  entitled  this  report  "Damping  of 
Very  Long  Internal  Gravity  Waves  in  the  Atmosphere"  (partly  carried  over 
from  F.Y.  1970). 


THE  W.K.B.  APPROXIMATION,  TURNING  POINTS  AND 
THE  MEASUREMENT  OP  PHASE  VELOCITIES 


1 .  Introduction 

In  the  fields  of  acoustic^*^  and  electromagnetic^*^wave  theory, 
as  well  as  in  some  aspects  of  elastic  wave  theory^,  K.K.li.  type 
methods  often  ^pfer  the  only  sensible  means  for  obtaining  useful, 
closed- form  approximations  for  the  description  of  wavefields  in  strati¬ 
fied  media.  The  method  is,  in  many  respects,  ideally  suited  for ’the 
calculation  of  the  far  field  in  wave  ducts^,  of  diffracted  fields 
under  conditions  of  caustic  formation1.  *  1 ,  of  focalization  intensities, 
etc.  Errors  incurred  by  the  use  of  this  method  arc  often  no  greater 
than  experimental  uncertainties,  such  as  those  due  to  fluctuations  of 
the  medium^;  indeed,  were  it  not  for  the  fact  that  the  approximation 
diverges  near  turning  points,  the  method  would  be  quite  sufficient  to 
deal  with  most  practical  problems  and  the  calculation  of  wavefields  by 
means  of  exact  solutions  of  the  wave  equation  would  acquire  a  largely 
eclectic  interest.  Thus,  in  long  range  transmission  calculations  in 
geoacoustics  and  in  e.m.  theory,  the  cumulative  experimental  uncertain¬ 
ties  of  phase  are  often  comparable  to  or  greater  than  the  errors  of  the 
N.k.B.  procedure:  as  long  as  one  can  stay  away  from  turning  point 
effects,  there  is  then  no  reason  to  try  to  improve  upon  this  approximation. 

There  exists,  however,  another  class  of  problems  in  which  one  tries 
to  make  accurate  phase  comparisons  over  short  ranges.  This  occurs,  for 
example,  in  the  design  or  interpretation  of  array  measurements  (i.c.,  of 
any  system  consisting  of  two  or  more  sensors  comparing  phase).  The  usual 
methods  of  relating  phases  for  neighboring  measurements  arc  equivalent  to 
the  tacit  assumption  that  the  medium  is  locally  homogeneous  or  that  ray 
optics  and  the  W. K.B.  picture  are  valid;  while  this  assumption  is  fre- 
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quvnth  justified,  it  can  and  will  break  down  in  stratified  geophysical 
media  in  the  neighborhood  of  turning  points.  Thus,  in  a  vertically 
stratified  medium  (in  which  the  physical  parameters  are  functions  of  the 
vertical  ;  coordinate),  this  breakdown  may  lead  to  false  conclusions 
concerning  vertical  phase  velocity  and  directionality.  These  observa- 
ti  ns,  which  we  shall  clarify  in  sections  3,  4  and  5,  are  germane,  in 
principle,  to  geoacoustic  fields  like  underwater  sound,  atmospheric 

i 

mfrasonics  and  c.m.  waves.  As  we  shall  see,  they  are  probably  mest 
pronounced  for  long  period  acoustic-gravity  waves  in  the  atmosphere. 


The  W.K.B.  method  and  the  turning  point  phenomenon:  review 
Consider  a  stratified  medium  in  which  the  physical  parameters  are 
functions  of  z  only.  Whether  we  are  dealing  with  acoustic,  gravity  or 
e.m.  waves,  the  assumption  of  simple  harmonic  wavetrains,  i.e.,  of 
solutions  proportional  to  e^°x  reduces  the  problem  to  a  Helmholtz 
equation  of  the  type: 
d‘ h 


dTT  4  Yfch  ■  0 
where  y  is  some  function  of  z. 

The  classic  W.K.B.  approximation  has  the  form: 

h  =  Y'Vi(s’So) 

with 

s  =  *  (  ydz 


CD 


(2) 


(3) 


and  s  is  an  integration  constant. 


As  is  well  known 
long  as: 


(1,2) 


,  equation  2  is  a  good  approximation  as 


1  ]  1  d2  .  2  1  [  d  . 

ar? lnY  •  lni,J 


«  1 


(4) 


s 


a  condition  that  is  seen  to  break  down  near 

Y  =  0  (5) 

i ,e. ,  at  a  turning  point  of  the  Helmholtz  equation  (1).  With  the 
exception  of  the  neighborhood  of  multiple  turning  points  (for  which  one 
or  more  derivatives  of  y  also  vanish),  (4)  is  equivalent  to: 

ij-£1„Y|  «  1  (6) 

This  somewhat  more  succinct  statement  says  that  the  W.K.B.  approximation 
applies  as  long  as  the  physical  properties  of  the  medium  do  not  change 
appreciably  over  a  vertical  wavelength;  it  may  also  be  shown^  that 
this  is  equivalent  to  assuming  that,  as  the  wave  energy  propagates 
through  the  medium,  one  may  neglect  the  process  of  continuous  backscattcr 
(reflection)  from  the  stratification.  Near  the  turning  point  itself,  the 
z  wavelength  becomes  very  large,  conditions  (4)  and  (6)  break  down  and  the 
approximation  (2)  becomes  divergent.  Exact  solutions  near  the  turning 
point  of  course  show  no  unusual  behavior  at  y  =  0:  this  point  simply 
marks  the  transition  from  a  region  of  oscillating  solutions  (y2  >  0)  to 
exponential  ones  (y2  <  0):  the  wave  function  and  all  its  derivatives  are 
continuous  through  the  turning  point. 

It  is  useful  and  instructive  to  examine  the  case  of  a  linear  y2: 

y2  *  Pz  ♦  y2  (7) 

in  which  case, 

s  ■  *3!  Y3  .  for  Y2  >  0  (8) 

V 

Standing  wave  or  modal  solutions  may  then  be  expressed  in  terms  of 
Bessel  functions: 

h  «  y£AJ  (s)  ♦  BJ  (s)]  •  (9) 

n  *6 

i-  (s>  • 


(10) 
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and,  in  the  y2  <  0  region: 

t  ?  -  2 

V  *-  —  "  v*- 


o  =  ‘Jit  Y 


•  3 


h  =  Y '  [-AI  (o)  ♦  BI  (a)] 


1/3 


-1/3 


(11) 


(12) 


77  =  Y*  2  [AI  (o)  -  BI  Co)) 
U-  -2/3  / 3 


(13) 


where  I  ,  I  „  are  the  modified  Bessel  functions. 

*  1 / 3  *2/3 

It  is  a  simple  matter  to  verify  that  the  W.K.B.  approximation  cor¬ 
responds  to  the  asymptotic  forms  of  (9),  (12)  for  large  s  or  a,  an 
observation  that  provides  a  direct  and  classic  method  for  determining 

the  connection  formulas  between  W.K.B.  solutions  to  both  sides  of  the 
(8) 

turning  point'-  ' . 

The  progressive  wave  solutions  corresponding  to  (9) ,  (10)  have  the 

form: 

( 1  91 

(14) 


h  *  Y»[J,2)(s) 

*/  3 


§=y2  H_l:r'c») 


,(1.2) 

-2/3 


(15) 


where  Mankel  functions  of  the  first  or  second  kind  are  used,  depending 
upon  whether  one  wishes  to  represent  down-  or  up-traveling  disturbances 
(depending  upon  the  sign  of  the  exponent  in  e*iu*).  This  is  conveniently 
seen  by  going  to  the  asymptotic  forms: 


lljl)  «  2is(ns)'V(s  "  5*/12) 

II *  2^( irs) ’  5n/12) 


(16) 

(17) 


he  see  that  the  W.K.B.  approximation  (2)  corresponds  to  the  use  of 
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those  asymptotic  approximations.  Liquations  (16),  (17)  are  in  fact  quite 
accurate  for  s  >>  1  -  a  condition  which,  by  virtue  of  (8)  and  (7)  is 
equivalent  to  the  statement  —  lny  <<  — ,  i.e.,  is  consistent  with  (6). 

In  dealing  with  a  medium  of  indefinite  extent  and  linear  y2,  we 
require  that  in  the  region  y2  <  0  solutions  (12),  (13)  give  a  convergent 
energy  density  for  a  ■*  <*>.  This  is  achieved  by  taking  A  =  B,  giving  for 
h: 

f 

h«A3\-V  K.  (o)  (18) 

73 

—  0 

where  k,  .  (o)  vanishes  like  o  e  for  large  o. 

1  '  3 

Note  that  the  standing  wave  solutions  (9)  with  A  =  B,  together  with 

(15) ,  are  the  result  of  superposing  an  incident  and  a  totally  reflected 
field.  Thus,  taking  the  turning  point  y  =  0  at  the  origin  z  =  0,  with 
y2  >  0  for  z  <  0,  we  have  P  <  0,  and  taking  s  *  -/z  ydz  =  -^j>73  gives 

s  >  0  for  y2  >0;  with  these  conventions,  the  incident  wavetrain  traveling 
towards  the  turning  point  would  be 

*i  =  A  ySl[V  (s)  g-iwt  (19) 

and  the  reflected  wave  train,  traveling  away  from  the  turning  point: 

*R  =  A  e1*  yHl^fs)  e’lwt  (20) 

It  is  easily  verified  that  the  superposition  of  (19)  and  (20)  gives  the 
standing  wavcfield  (9)  with  A  =  B,  and  gives  (18)  in  the  y2  <  0  region. 

It  is  important  to  understand  here  that  this  total  reflection  process 
is  spread  over  a  zone  near  y  =  0.  Thus,  as  s  ■*  0  we  approach  the  turning 

iJL 

point  and  ■+  -c  <K  ,  so  that  the  actual  reflection  coefficient  at 

H  -2i4  * 

z  =  0  is  equal  to  -e  =  e  ,  on  the  other  hand,  as  s  -*■  °°,  wc  see  from 

(16) ,  (17)  that  -*■  e^^s  i.e.,  the  effective  reflection  coef¬ 

ficient  is  asymptotically  equal  to  e"^.  There  is  then  a  region,  in  the 


VI.  I  II  I  l  V  .'I 


tin-  turning  point,  tli.it  t  .ml  ri  hull’*,  uu  ovtin  phase  change  o  J 


*  o,  above  and  beyond  the  integrated  phase  path  contribution:  in  other 
words,  the  reflection  process  is  not  concentrated  at  the  turning  point, 
but  is  smeared  out  over  a  layer  of  thickness  given  by  s  -  1. 

It  is  also  worth  noting  that  while  h  and  dh/dz  both  change  from 


oscillating  to  exponential  type  solutions,  they  do  not  obey  the  same 
C9i 

Helmholtz  equation  .  Thus,  differentiating  (1)  gives: 


,dh  ,  dy2 

♦  y2T“  +  h-f— 
'  dz  dz 


Taking 


It  appears  then  that  W  satisfies  the  Helmholtz  equation  (25)  with  a 
turning  point  at 

6  -  0  (27) 

i.e.,  at  a  different  value  of  z  from  that  corresponding  to  (s) . 

One  might  conclude  from  this  that  the  function  y  =  dh/dz,  or  some 
linear  combination  of  h  and  y,  changes  character  (from  oscillating  to 
exponential)  at  a  different  value  of  z  than  h  itself.  That  such  a  con- 
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elusion  i s (  in  fact,  quite  specious  is  seen  from  the  exact  solutions 
(y) ,  (10)  and  (12),  (13).  Mathematically,  tills  apparent  ambiguity  arises 
from  the  feature  that,  at  y  =  0,  R2  is  singular:  a  proper  discussion  of 
(25)  near  the  singularity  must,  needless  to  say,  yield  the  same  conclu¬ 
sion  as  examination  of  the  exact  solutions. 

It  follows  also  from  what  has  been  said  above  that  the  ray  picture 
interpretation  of  y,  B  is  misleading  near  turning  points.  Conclusions 

i 

concerning  direction  of  phase  travel,  energy  flux,  levels  of  reflection, 
etc.,  that  are  based  upon  one's  usual  intuitive  picture  of  the  wave- 
number  vector  arc,  in  the  final  analysis,  only  as  good  as  the  K.K.B. 
approximation;  furthermore,  inspection  of  (4)  and  (26)  shows  that  the 
domain  of  validity  of  t hi s  approximation  is  precisely  the  same  for  h,  y, 
w,  dw/d:  or  for  any  linear  combination  of  these  functions. 

We  shall  show  in  what  follows  that  there  are  situations  in  geo- 
acoustics  and,  in  particular,  in  acoustic-gravity  wave  theory,  for  which 
these  features  of  the  wave  solutions  acquire  a  practical  significance. 

5.  Phase  velocity  measurements  in  stratified  media 

We  now  focus  our  attention  on  the  problem  of  measuring  the  z  compon¬ 
ent  of  phase  velocity  (normal  to  the  stratification),  a  problem  one 
encounters  in  the  use  of  vertical  arrays  in  the  oceans  or  atmosphere. 

This  situation  occurs,  for  example,  whenever  one  attempts  to  interpret 
vertical  phase  velocities  in  terms  of  the  stratification  (structure)  of 
the  medium. 

A  general  representation  of  a  traveling  harmonic  wavefield  is: 
i  d>  -  i  uit 

$  =  m  e  v  e  .  (28) 

where  m,i£  are  real  functions.  The  phase  velocity  in  the  z  direction  is: 

0) 

v  a  - - - 

z  3$/8z 


(29) 
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1 

\ 


It  (2)  is  valid,  it  follows  from  (3)  that: 

3^/Dz  =  ±Y 
and 


to 

v  =  *— 

Z  Y 


(30) 


(31) 


Thus,  if  the  W.K.B.  approximation  applies  our  usual  understanding  of 
vertical  phase  velocity  as  the  ratio  of  angular  frequency  to  vertical 
wavenumber  component  is  valid:  this  is  the  realm  of  ray  optics.  , 

But  if  the  conditions  of  validity  for  the  W.K.B.  method,  i.e., 
inequalities  (4)  or  (6),  are  not  fulfilled,  the  statement  (31)  is  no 
longer  correct.  This  is  easily  seen  by  turning  one's  attention  once 
again  to  the  exact  solutions  for  the  linear  Y ~  case.  Here  $  is  given  by 
the  phase  of  H^'^(s),  i.e., 


!/ 3 


$  =  arctan  \  cotan^- 


i_.  J->A1(*)  1 


*  7T  J,,  (S 

smj  */  3 v 


(32) 


giving 


li  =  t  il 

9z  tts 


Ji/3(s)  +  i[Jv3(s)  -  2J-i/3(s)]2^ 


(33) 


By  going  to  the  limit  of  vanishing  s,  y  (using  equation  8),  we  have  the 
vertical  phase  velocity  at  the  turning  point: 

as  Y  -  °.  vz  =  ->  t4n  3'V3|P|‘1/3[r(2/3)]-2UJ  (34) 

Thus  the  vertical  phase  velocity  tends  to  a  finite  limit  at  the  turning 

point  (whereas  literal  acceptance  of  the  ray  picture  leads  to  an  infinite 

value  for  v  ) . 

zJ 

Therefore,  unless  exact  solutions  for  41,  9 4>/ 9z  are  used,  one  may  be 
appreciably  in  error  in  deducing  the  vertical  phase  velocity  and  angle  of 
travel  of  a  harmonic  wavefront.  Conversely,  for  measurements  made  near  a 
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turning  point  interpretation  of  observed  phase  velocities  may  lead  to  er¬ 
roneous  conclusions  cc  rning  the  stratification  unless  the  exact  equa¬ 
tion  (33)  is  used  in  ). 

Figure  1  shows  a  it  of  $  and  of  d4>/ds  =  •d'J'/dz.  These  quanti¬ 

ties  are  plotted  aga  s,  so  that  the  W.K.B.  oximation  for  <p  is, 
by  (2),  a  straigh*  of  slope  one;  it's  'pt  is,  by  (16),  s0  =  5tt/12. 

The  exact  curve  o 4>/ 3 s  are  seen  to  a  ''he  W.K.B.  val.ues 

quite  closely  £oj 

s  »  1  (35) 

At  s  =  1,  the  W.  1  rror  for  3<fc/3s  and,  thu.;.,  tical  phase  velo¬ 
city  calculation  less  than  5?»:  condition  (i  as  pointed  out 

above,  an  altem-  statement  of  (6). 

It  is  <  ue  interest  to  estimate  (34)  low  frequency 

acoustic  waves  (v  ill  discuss  acoustic-gravity  ?  in  the  next 

section).  If  on  lects  the  effect  of  gravity  v  substitute  for 

P  in  (34)(1): 

P  =  oj  ^  (36) 

where  c  is  the  s<  peed. 

This  gives  t  rning  point  value: 

\  a  2  1/3f1/3  (37) 

c 

where  f  is  the  f  zy  in  hz.  In  underwater  acou'  one  usually 

.  1/, 

deals  with  gradients  less  than  dc/dz  'v  1  sec' 1 .  Thus  2.3  f  3  is  pretty 
much  the  smallest  value  for  (37)  in  marine  acoustics  experiments.  It 
gives  v  /c  -  5  for  f  =  10  hz,  and  v  /c  =  10  for  f  =  102  hz.  In  practical 
terms,  the  difference  between  an  infinite  move-out  velocity  and  one  that 
takes  place  at  10  times  the  speed  of  sound  is  hardly  significant,  when 
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oiu  i\  iombcrs  that  the  mcusui'cuicnt  must  take  plucc  over  a  fraction  of  a 
wa\e  length;  in  underwater  sound  experiments,  for  frequencies  of  102  hz 
and  above,  this  effect  is  unlikely  to  be  important. 

However,  in  infrasonic  measurements  in  the  atmosphere  and  thermo¬ 
sphere,  the  phenomenon  is  more  likely  to  be  observable.  Thus,  gradients 
of  the  order  of  dc/d:  =  10-2  sec-1  are  characteristic  of  the  thermosphere. 
Signals  having  periods  as  long  as  4  minutes  have  been  observed  in  con- 
nectioa  with  Apollo  firings  (10),  i.e.,  we  may  take  f  2  4*10‘3  and 
v../c  2  at  the  turning  point:  here  measurements  of  the  critical  phase 

velocity  may  yield  misleading  results,  unless  they  are  interpreted  with 
the  help  of  exact  formulas  of  the  type  of  (29),  (33).  As  we  shall  see  in 
Section  4,  this  effect  is  even  more  marked  for  acoustic-gravity  waves  in 
the  5-30  minute  period  bandpass  at  thermospheric  heights. 

As  for  the  behavior  of  dh/dz  near  the  turning  point,,  we  once  again 
have  recourse  to  the  exact  form  for  the  progressive  wave  solution  with 
linear  y2,  i.e.,  we  use  equation  15.  Thus  dh/d*.  corresponds  to  wavefields 
of  the  type: 


...  i^  -iwt 

V  =  n  e  c 


where  the  phase  1;  is,  by  (20) 


TT  J  ^2/3 

0  =  arctan  [ cotan  •=•  ♦  — : - r=-  •  ~ - 7-r] 

1  3  sin  ir/3  J.2/3^ 


(38) 

(39) 


For  large  s,  i.e.,  in  the  region  for  which  W.K.B.  applies,  it  is  clear 

f  1  2) 

from  the  usual  asymptotic  form  for  (s),  that  ^  -*■  s  *  tt/  12  so  that 

jj/hz  -  3$/3z  and  the  phase  velocities  of  the  4>,  V  fields  are  the  same. 
But,  for  s  <  1,  equation  (39)  shows  that  the  behavior  of  3^/3z,  3$/3z 
will  be  quite  different.  In  fact,  if  we  use  the,  first  terms  of  the 
power  series  for  J+2/3(s)»  equation  (39)  gives 

1/3 


^  =0.15  s 
ds 


(40) 
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thus  the  phase  velocity  near  the  turning  point  goes  to  infinity  like  'Cl. 
The  behavior  of  v*  dv/ds  is  shown  in  I-'iguro  1. 

Acoustic-gravity  and  internal  gravity  waves 

he  have  established  that  one's  interpretation  of  infrasonic  wave 
behavior  near  turning  points  may  be  seriously  in  error,  unless  he  uses 
the  exact  field  equations.  The  use  of  ray-optical  concepts  here  is 
erroneous;  furthermore,  the  amount  and  sign  of  the  error  depends  upon 
the  field  quantity  being  studied. 

The  largest  errors  will  be  encountered  for  acoustic-gravity  waves 
at  ionospheric  heights.  Measurements  of  traveling  wave  systems  at  those 
heights  are  currently  a  subject  of  considerable  interest;  naturally 
occuring  wave  systems,  in  particular,  are  being  studied  intensively  for 
the  light  they  shed  upon  upper  atmospheric  structure^11*  12^;  disper¬ 
sive  acoustic-gravity  waves  with  periods  of  several  minutes  have  been 

detected  by  vertical  phase  sounders  in  connection  with  Apollo  launches 

fl3) 

and  disturbances  in  the  same  passbandv  '  as  well  as  with  longer 
periods  have  been  seen  in  connection  with  large  atmospheric 

explosions.  It  is  therefore  important  to  sec  how  the  considerations  of 
Section  3  apply  to  acoustic-gravity  and  internal  gravity  waves  at  these 
heights . 

V.'ith  this  in  mind,  let  us  fix  our  attention  upon  a  specific  model 
of  a  specific  section  of  the  thermosphere.  Since  much  of  the  recent 
phase  sounder  work  has  been  conducted  with  sounding  frequencies  in  the 
2-o  Mil z  passband,  i.c.,  at  true  daytime  heights  in  the  ISO  -  250  km 
range,  let  us  select  a  reference  height  of  180  km  and  use  this  as  the 
origin  of  our  z  coordinate.  We  may  then  expand  y2  into  a  Taylor  scries 
in  z:  keeping  the  first  term  only  gives  the  linear  law  (7).  The  range 
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o:  validity  of  this  model  depends  upon  the  actual  rates  of  variation  of 
the  physical  parameters.  We  begin  with  tho  differential  equation  for 
the  elastic  pressure  p^“^.  Writing: 


p  3  P 


1/2  h  ei(«x  -  ut) 


gives  equation  1,  with 


y2  3  -  u2  )  *  — r(.\’2  -  u>2) 

C*-  op  v~ 


(41) 


(42) 


where  c  is  the  adiabatic  sound  speed,  0  the  density,  u  is  the  elastic 
pressure  "acoustic  cutoff",  X  is  the  Brunt-VaislilU  frequency,  the  hori¬ 
zontal  phase  velocity 

It  can  be  shown 


v  ~  — 
x  a 


(43) 


that  a  reasonable  fit  for  the  actual  u  ,  c, 

op*  * 

N  variations  in  the  thermosphere  is  obtained  by  taking  the  perfect  gas 

law  with 


:2  =  Pz 


♦  c 


i.O.  , 


N"  =  ^7 (PZ  ♦  0.4g) 


1 


,2t2 


P 2) 


(44) 

(45) 

(46) 


"op  = 

where  t  =  1.4  is  the  ratio  of  specific  heats.  We  may  therefore  approxi¬ 


mate  y  by  a  linear  law  as  long  as  1/c2  is  approximately  linear,  i.e.. 


j : 

4  .  Ad  - 

c2  C  *  C 

0  0 

The  ex.»c t  solutions  of  Sections  2  and  3  will  thus  be  useful  if 

■-=-  1 

c 

o 

Taking  at  :  =  0  (ISO  km  height)  the  values^). 

c  2  =  5.94*105  m2  sec"2 
0 

P  -  6.2  m  sec"2 


(47) 


(48) 


(49) 
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*  *  l  m  tl;a ;  (48)  imposes  the  limits: 

j=!  <<  ID5  a  (50) 

la  other  words,  we  shall  be  limited  to  a  range  of  about  *25  km  to  each 
si«.o  of  -  =  0,  i.e.,  to  a  zone  of  true  heights  between  155  km  and  205  kn. 
Kith  this  limitation  in  mind,  we  can  apply  the  results  of  Sections  2  and 
3  to  the  actual  thermosphere, 
liquations  (45),  (47),  (49)  give  then: 

X2  =  Nq2  ♦  n:,  n  =  -  *^02“7’  (51) 


Using  the  numerical  values  (49)  and  (51),  (52)  we  may  calculate  exact 
wave  functions,  phases,  phase  velocities  etc.,  for  155  km  <  z  <  205  km. 

The  procedure  may  be  applied  in  the  vicinity  of  a  succession  of  z  values 
to  cover  any  desired  height  range  in  the  thermosphere. 

Figures  2  and  3  show  the  behavior  of  s  at  these  heights,  for  the 
thermospheric  constants  (49),  (51),  (52):  it  is  seen  that  the  validity  of 
the  W.K.B.  and  ray-optical  points  of  view  is  contingent  upon  frequency  and 
v  . 

X 

The  vertical  displacement  (  is,  to  within  a  ±  iu  factor,  similar  to 
the  vertical  velocity  and  is  of  particular  interest  since  it  is  more  direct¬ 
ly  related  to  the  physical  quantity  measured  in  ionospheric  phase  sounders 


* 


l*.*  Sod  r.Hi  velocity),  It  is  related  to  the  elastic  pressure  l'ield^1*^  by 
the  equation: 

c  •  (^i5  *  i&) 1  p(u2  •  &  (55) 


But 


3n  Id  .  i/2  dh  i(ax  -  ut) 

'  P  T  XT  lnp  ♦  0  '  nT  c 


J;  r  2  il: 


(56) 


and 


«2.-l 


t  =  p’l/2(u2  -  {=7)’  [  h(i  35.  lnp  ♦  |r)  ♦  jj) 


v' 

X 


(57) 


’ll  10  relative  importance  of  the  h,  hz  terms  near  the  turning  point  is 
determined  by  examining  the  exact  solutions  (14),  (15)  in  the  limit  of 
Y  -•  0,  s  ■»  0.  It  is  then  seen  that  dli/dz  a  0.5  P^3h.  But,  for  the  case 
under  consideration,  for  periods  in  the  10  -  30  minute  range  and  for  hori¬ 
zontal  phase  velocities  vx  <  300  m  sec"1.  P  will  be  seen  to  vary  in  the 


10 


- ;  3 


to  10*i,4m"J  range.  Sinco  furthermore  %/v*  >  10’4  m”1  and 


,-iV3 

1"  5  m-l 


d/d;  lno  ~  10  3  m  we  have,  approximately: 


;  =  0-1/2  («r  -  -4)-‘  Ah  (58) 

X  X 

This  approximation,  good  for  v^  <  3X102  m  sec’1  and  periods  in  the  10  - 
30  minute  range,  shows  that  we  may  usefully  apply  the  exact  solutions  (14) 
to  obtain  the  vertical  phase  velocity  from  (29)  and  (33)  and,  in  particular, 
the  limiting  turning  point  value  (34).  Since  we  have  here  vx  /c  <<:  1,  it 
follows  from  (53)  that  we  have,  approximately: 

N  2  P 

(59) 

X  o 

and  the  turning  point  phase  velocity  (34)  is  given  by 


s  4tt  3-4/3  [r(4)]’2  N-2/3  P-1/3c2/3wv-1/3 
v  1  v3  J  ox 

x 

a  1 ,4*10J*tuxv'1/3 
X 


(60) 
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Thu-;,  for  periods  of  the  order  of  30  minutes  and  horizontal  phase  velocities 
~  Id-  n  sec”1  we  have  v,/vx  ~  1  at  the  turning  point,  giving  an  apparent 
phase  front  inclination  of  45° !  Since  v,  -  vx  is 

often  observed  in  the  ionosphere ,  these  calculations  are  very  perti¬ 
nent  to  the  interpretation  of  actual  ionospheric  observations  of  acoustic 
and  acoustic- gravity  travelling  wavetrains. 

It  must  be  emphasized,  however,  that  in  general  it  is  (57)  that  must 
be  used:  the  interpretation  of  phase  velocity  measurements  must  then 
proceed  by  defining,  with  the  help  of  (14)  and  (15),  an  exact  phase  function 
for  h.  This  is  slightly  more  involved  than  the  calculations  for  our  sim¬ 
plified  case  with  vx  <  3*102  m  sec"1,  but  does  not  offer  any  particular 
di fficulties. 


5.  Discussion 

In  problems  involving  wave  propagation  in  stratified  media,  it  is  not 

difficult  to  convert  a  set  of  values  for  the  local  parameters  and  their 

derivatives  into  an  explicit  ,  simple  criterion  for  the  validity  of  the 

K.K.B.  approximation  in  the  solution  of  the  fundamental  Helmholtz  equation 

(1) .  It  is  rather  illuminating,  in  fact,  to  assume  that  Y2  in  equation  (1) 

may  be  treated  locally  as  a  linear  function  of  z  (a  classic  assumption 

(8) 

first  used  by  Langer'-  •'in  his  solution  of  the  connection  problem,  the  limi¬ 
tations  of  which  are  easily  stated  in  terms  of  a  permissible  range  of  z 
values).  This  approach  has  the  advantage  of  yielding  locally  exact  solu¬ 
tions  of  (1)  that  are  valid  for  simple  turning  points  y  =  0,  solutions 
which  furthermore  yield  the  W.K.B.  approximations  explicitly  as  an  asymp¬ 
totic  limit.  This  shows  that  the  W.K.B.  method  is  valid  in  the  limit  of 
large  enough  phase  paths,  a  statement  that  can  be  written  very  succinctly 
as  s  >>  1,  where  s  is  given  by  (3).  The  W.K.B.  approximation  is  essentially 
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the  piano  wavo  ray-optical  point  of  view,  so  that  this  criterion  yields 
also  the  domain  of  validity  for  the  usual  formulas  relating  phas<  velocity 
(vertical  phase  volocity,  in  the  present  case)  to  the  corresponding  wave¬ 
number  components  and  angles  of  incidence.  In  fact,  a  glance  at  Figure  1 
shows  that  the  criterion  may  bo,  in  practice,  rolaxcd  to  s  >  1;  ;n  this 
domain  everything  can  be  interpreted  in  torms  of  the  usual  ray-optical 
imagery.  In  the  vicinity  of  a  turning  point  of  the  Helmholtz  equation, 
interpretation  becomes  more  delicate:  the  W.K.B.  approximation  is  no  longer 
valid  and  the  usual  concepts  of  wavenumbor  and  phase  velocity  components 
tail  us  -  but  wo  do,  in  this  case,  have  exact  equations  to  fall  back  on. 

The  failure  of  the  simple-mindod  ray-optical  concepts  here  is  best  under¬ 
stood  as  the  result  of  interference  effocts  concentrated  near  the  turning 
point.  Thus,  we  may  visualize  the  flow  of  wave  energy  through  a  continu¬ 
ously  stratifiod  medium  as  actually  taking  place  through  a  succession  of 
very  thin  homogeneous  layors,  representing  incremental  changes  in  the  physi¬ 
cal  parameters.  Assume,  for  the  purpose  of  simple  visualization,  that  we 
arc  dealing  with  acoustic  waves  in  a  temperature  stratified  medium;  when 
the  wavefronts  and  wavenumber  vectors  can  bo  unambiguously  defined  they  are 
mutually  orthogonal,  and  the  energy  flux  vector  is  colincar  with  the  wave- 
number:  as  the  wave  travels  towards  regions  of  higher  sound  velocity,  the 
rays  (energy  paths)  are  bent  according  to  Snell's  law;  as  the  wave  progresses 
reflections  take  place  at  the  boundary  of  each  incremental  layer:  as  the 
layers  become  infinitesimal,  this  provides  a  continous  reflection  process, 
a  sort  of  effective  backscattcr  of  tho  energy.  As  has  been  pointed  out^, 
the  W.K.B.  approximation  neglects  this  backscattcr  entirely,  a  procedure 
that  is  justified  as  long  as  the  parameters  of  the  medium  vary  little  over 
a  vertical  wavelength.  But  then,  as  we  approach  a  turning  point,  we  near 
a  condition  of  grazing  incidence  so  that,  no  matter  how  small  the  increment 
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in  sound  velocity  botween  neighboring  layers,  one  eventually  approximates 
a  condition  of  total  reflection:  the  backscattcr  is  then  no  longer  small 
and  cannot  be  neglected.  This  is  the  physical  reason  for  which  the  W.K.B. 
method  fails;  it  is  also  the  reason  for  which  it  becomes  impossible  to  de¬ 
fine  a  vertical  phase  velocity:  even  an  isolated  "upward"  or  "downward" 
moving  wavefield  of  the  type  (19)  or  (20)  will  thus  get  involved  with 
interference  effects  near  y  =  0.  These  effects  are  irreducible  and  cannot 
be  eliminated  by  directional  filtering.  It  must  be  emphasized  that  this 
is  not  the  same  thing  as  the  generation  of  a  standing  wave  pattern  such  as 
(9)  (e.g.,  with  A  =  B)  by  a  superposition  of  the  incident  and  reflected 
wavefields  (19),  (20),  a  pattern  that  is,  in  principle,  analyzable  into  two 
separate  upward  and  downward  moving  wavefields.  It  is,  for  instance,  pos¬ 
sible  to  produce  a  pure  wavefield  of  the  type  (s)  e’*Wt  moving  away 

from  the  turning  point  (o.g.,  by  leakage  of  a  harmonic  wave  system  through 
a  high  velocity  barrier);  such  a  wave  would  be  carrying  energy  away  from 
the  turning  point  but  would,  nevertheless,  scatter  and  interfere  with  it¬ 
self  in  a  narrow  region  near  y  =  0.  Since  this  scatter  actually  consists 

(2  7) 

of  multiple  reflections  *  ,  it  contains  both  upward  and  downward  moving 

components,  so  that  directional  filtering  cannot  get  rid  of  these  effects 
which  are  thus  irreducible . 

This,  then,  is  the  root  of  the  difficulties  one  may  have  in  interpreting 
the  phase  velocities  of  a  travelling  harmonic  wavefield  in  the  vicinity  of 
turning  points  in  stratified  media,  i.e.,  in  regions  for  which  the  W.K.B. 
method  fails.  We  have,  in  Sections  3  and  4,  dealt  with  these  questions 
explicitly  in  the  case  of  infrasonic  and  acoustic-gravity  wave  propagation 
in  the  thermosphere.  In  the  last  few  years,  observations  have  been  made  of 
travelling  wave  systems  at  these  heights  and,  in  some  cases  the  ob¬ 

servations  appear  to  fall  in  the  range  of  values  for  which  interpretation 
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becomes  ambiguous.  Tho  purpose  of  this  paper  has  been,  chiefly,  to  explore 
and  clarify  tho  naturo  of  tho  ambiguity  in  such  a  way  that  future  experi¬ 
mental  measurements  at  inf "asonic  and  acoustic-gravity  wave  frequencies 
will  be  able  to  deal  with  it;  since  this  is  largely  a  matter  of  making 
sure  that,  in  such  instance,  one  has  enough  vertical  measurement  points 
to  determine  whether  or  not  he  is  approaching  a  turning  point  region,  this 
should  not  be  difficult  to  achieve.  , 
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Figure  Captions 

Figure  1.  Progressive  wave  phase  angles  as  a  function  s;  $  is  the  phase 
of  h,  4'  the  phase  of  dh/dz  in  a  medium  of  linear  y2(z)«  The  asymptotes  of 
4> »  v  correspond  to  the  W.K.B.  approximation.  It  is  seen  that  the  slopes 
<l>'  have  quite  different  behavior  near  s  =  0  (turning  point),  implying 
different  limiting  vertical  phase  velocities.  It  is  also  seen  that  phase 
and  phase  velocity  are  satis factorily  given  by  W.K.B.  for  s  >  1. 

Figure  2.  Behavior  of  s  in  tlv  earth's  thermosphere  at  an  altitude  of  180 
km,  as  a  function  of  u  for  various  selected  values  of  the  horizontal  phase 
velocity  v^.  The  dotted  line  s  =  1  separates,  for  all  practical  purposes, 
the  region  in  which  the  W.K.B.  approximation  is  valid  (s  >  1)  from  the 
region  in  which  it  is  insufficiently  accurate  for  the  calculation  of  phase 
velocity  (s  <  1);  note  that,  in  the  latter  region,  it  may  still  be  possible 
to  use  W.K.B.  for  order  of  magnitude  amplitude  calculations. 

Figure  3.  Same  as  Figure  2,  but  with  horizontal  phase  velocity  v^  as 
variable  and  period  in  minutes  as  parameter. 
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DAMPING  0I:  VliRY  LONG  INTLKNAL  GRAVITY 
KAVliS  IN  Tllii  ATMOSI'I  UiRl: 


Classic  waveform  calculations  for  utmosphoric  explosions,  such  as 
those  of  iiarkridcr  (1904),  have  been  largely  (but  not  entirely,  sec  c.g., 
Iiarkridcr  and  Wells,  1968)  limited  to  the  acoustic  and  shorter  period 
internal  gravity  wave  spoctrum  (periods  T  <  5  min.).  These  calculations 
reproduce  quite  well  the  waveforms  observed  on  microbarographs.  There 
is  now  a  growing  body  of  evidence  indicating  that  long  wavelength  distur¬ 
bances  generated  by  nuclear  explosions  in  the  atmosphere  can  also  travel 
to  great  distances.  Both  ionospheric  measurements  (Breitling  ct  al, 

1967;  Rose  ct  al,  1961;  Hultquist  et  al,  1961;  Dicmingcr  ct  al,  1962; 

Herron  and  Montes,  1970)  and  long  period  barographs  (Tolstoy  and  lierron, 
1970;  Herron  and  Montes,  1970)  indicate  the  occasional  existence  of  long 
wavelength  arrivals  having  travelled  at  least  once  around  the  globe. 
Propagation  to  these  great  distances  implies  that,  for  long  wavelengths, 
attenuation  need  not  be  prohibitive  and  that  channelling  and  waveguide 
effects  must  take  place.  Although  it  is  possible  that  the  surface  gravity 
node  (m  =  6)  may  account  for  some  of  these  observations  (Tolstoy  and 
Herron,  1970;  Tolstoy  and  Pan,  1970),  multiply  reflected  internal  gravity 
waves  arc  potentially  capable  of  explaining  arrivals  with  group  velocities 
below  600  m.sec-1.  Thus  we  shall  show  in  what  follows  that  the  attenuation 
of  the  lower  modes  (m  ■  1,  2,  3)  need  not  be  excessive.  The  relatively 
short  range  observations  of  Stoffrcgcn  ct  al  (1961)  can  be  explained  by 
internal  gravity  waves  (Hines,  1967)  along  direct  ray  paths  only, traveling 
from  the  explosion  source  to  the  observation  point,  with  no  reflections. 

It  has  been  suggested  that  these  waves  are  absorbed  in  the  upper  atmosphere 
because  of  what  amounts  to  an  indefinite  increase  in  kinematic  viscosity 
with  height  as  the  density  tends  to  zero.  In  other  words,  upwards  traveling 
energy  is  attenuated  and  transformed  into  heat:  there  is  no  reflected  wave 

24 


i 


energy.  Hi  is  situation  would  correspond  schematically  to  that  of  figure  1. 
As  wo  shall  sec,  this  mechanism  is  only  plausible  for  sufficiently  short 
wavelengths:  long  waves  arc  reflected  in  the  manner  demonstrated  below, 
an  effect  which  was  first  shown  by  Ynnowitch  (1967).  This  appears  to 
be  the  only  reasonable  explanation  for  the  observed  propagation  of  long 
wavelengths  to  ranges  of  10**  km  or  more. 


In  order  to  keep  the  discussion  simple,  we  shall  consider  a  locally 
isothermal  atmosphere,  i.e.,  an  atmosphere  consisting  of  layers  in 
which: 

■2vz 


P  =  PG  e 


(1) 


where  v  is  a  constant.  At  ionospheric  heights  between  200  and  500  km.. 


v  is  probably  such  that: 


1  x  l(f5  nf1  <  v  <  3  x  10"5  in-1 


(2) 


Assuming  then 


v  *  1.5  x  10-5  m'1 


(3) 


implies  that  the  thickness  d  over  which  the  density  varies  by  one  order 
of  magnitude  (the  scale  height)  is: 


d  »  2.3/2v  *  .8  x  105  m 


(4) 


Essentially  complete  absorption  of  the  energy  should  take  place  for  wave¬ 
lengths  of  the  order  a  few  d  and  less:  this  covers,  in  practice,  a  good 
part  of  the  observed  TID  spectrum. 

The  statement  that  the  energy  of  an  upward  traveling  wave  gets 
entirely  absorbed  and  changed  into  heat  hinges  on  a  literal  acceptance 
of  the  equations  of  continuum  mechanics  and  upon  the  ensuing  infinite 
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kinematic  viscosity  as  p  ►  0.  This  statement  must  also  he  qualified  in 
terms  of  the  wavelength:  it  has  been  shown  by  Yanowitch  ( 1 067)  that  waves 
of  sufficient  length  arc  reflected;  Yanowitch  interprets  this  result  to 
mean  that  one  has  reflections  from  layers  of  large,  rapidly  varying  kine¬ 
matic  viscosity.  Yanowitch's  proof  is  based  upon  a  straightforward  dis¬ 
cussion  of  the  effects  of  the  viscous  terms  on  the  internal  gravity  wave 
field  equations. 

In  fact,  of  course,  the  continuum  equations  cease  to  be  valid  as 
o  *  0,  and  a  more  sophisticated  treatment  is  required:  needed  here  is  an 
adequate  molecular  model  of  the  gas,  combined  with  an  analysis  of  Boltz¬ 
mann's  equations  for  a  rarefied,  density  stratified  gas  in  a  gravity  field. 
Such  an  analysis  would  be  extremely  difficult  and  is  beyond  our  reach  at 
present,  however,  pending  the  appearance  of  a  more  rigorous  theory  along 
such  lines,  one  may  perhaps  indulge  in  some  assumptions  as  to  what  the 
long  wavelength  limit  of  such  solutions  is  likely  to  show  (Tolstoy,  1967). 
Interestingly  enough,  it  is  possible  in  this  manner  to  obtain  a  result 
comparable  to  Yanowitch's. 

\ 

We  assume  the  critical  dimensionless  parameter  in  this  problem  to  be 
kn,  where  k  is  the  wavenumber  and  i  is  the  mean  free  path  of  the  molecules. 
For  angles  that  arc  not  too  steep,  it  would  seem  plausible  to  assume  that: 

1)  For  heights  such  that 

kfc  <<  1  (5) 

the  equations  of  continuum  mechanics  will  apply 

2)  For  heights  such  that 

kt  >  1  (6) 

i.e. ,  if  the  mean  free  path  of  the  molecules  is  of  the  order  of  or  less 
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than  the  wavelength,  the  medium  acts  effectively  as  a  vacuum  [providing 
hydromagnet ic  interactions  can  be  neglected  in  the  atmosphere  at  these 
heights:  this  is  probably  so  for  periods  less  than  three  hours  or  so, 
as  can  be  shown  by  application  of  the  Dungcy-Fejcr-llines  criterion, 
see  e.g.,  Tolstoy,  1967].  An  equivalent  statement  would  be  that,  for 
collision  times  longer  than  the  characteristic  time  scale  (period)  of 
the  disturbance,  no  transmission  of  energy  can  occur  (this  ic,  effective 
ly,  the  operational  definition  of  a  vacuum). 

3)  Between  the  vacuum  region  defined  by  (6)  and  the  continuum  corres¬ 
ponding  to  (5)  there  exists  a  zone  of  high  kinematic  viscosity  in  which 
most  of  the  attenuation  takes  place.  We  shall  refer  to  this  as  the 
transition  layer.  The  height  of  this  layer  above  the  earth's  surface 
depends  upon  the  wavelengths  considered:  according  to  our  best  estimate 
of  the  mean  free  path,  the  base  of  the  layer  should  be  at  a  height  of 
the  order  of  500  i  50  km.  for  wavelengths  in  the  5xl02  -  5xl03  km.  band. 

Kinetic  gas  theory  tells  us  that  the  kinematic  viscosity  n  is  given 
approximately  by: 

n  =  3  <v>  i  (7) 

where  <v>  is  the  mean  velocity  of  the  molecules  and  l  is  their  mean  free 
path.  The  same  theory  tells  us  that 


<v>  =  c  (8) 

where  c  is  the  adiabatic  sound  speed. 

Elementary  wave  theory  further  shows  that,  for  constant  n,  a  pro- 
gressive  harmonic  wave  in  a  continuous  medium  attenuates  like  e  ^x,  where 


q  =  nk2/v 


ph 


(9) 


Vph  being  the  phase  velocity.  If,  for  an  order  of  magnitude  argument, 
we  assume  that  Vp^  and  c  are  of  the  same  order,  i.e. 
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V  »  c  =  <v> 

ph 

(10) 

and 

q  =  l  (k*)k 

(ID 

thus , 

for  fixed  kfc,  as  k  ■*  0,  q  -►  0. 

Let  us  then  assume  we  have  free  space  for 

k£,  >  1  (12) 

and  a  continuum  for 

* 

k i  <  10_1  (13) 

This  is  tantamount  to  the  statement  that,  for  given  k,  the  attenuation 
takes  place  primarily  within  a  layer  in  which  l  varies  by  one  order  of 
magnitude;  since  p  £  l,  it  follows  that  the  thickness  of  the  transi¬ 
tion  layer  is  essentially  the  scale  height  d  (equation  4). 

Finally,  we  make  the  additional  assumption  that,  within  this  tran¬ 
sition  layer,  the  attenuation  for  very  long  wavelengths  is  similar  to 
that  of  a  layer  of  thickness  d  having  a  mean  kinematic  viscosity  n: 

n  =  n°  ^  f  e2vz  dz  =  4n<>  (14) 

‘  o 

where  no  is  the  kinematic  viscosity  at  the  bottom  of  the  transition  layer. 

In  this  picture,  a  long  wavelength  (k  0)  internal  gravity  wave 
traveling  upwards  suffers  some  attenuation  as  it  passes  through  the  tran¬ 
sition  layer,  but  is  then  partially  reflected  (figure  2):  it  can  be 
plausibly  assumed  that,  for  long  enough  wavelengths,  this  reflection  pro¬ 
cess  is  similar  to  that  obtaining  at  a  "free"  surface  of  contact  between 
a  continuous  fluid  and  a  vacuum.  In  the  present  model,  interposed  between 
the  vacuum  and  the  continuous  fluid  is  the  transition  layer;  the  thickness 
of  this  layer  is  of  the  order  of  a  scale  height  and  is,  in  an  isothermal 
model  of  the  upper  atmosphere,  independent  of  the  wavelength:  for  very 
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long  wavelengths ,  therefore,  it  acts  effectively  as  a  thin  high  viscosity 
film. 

The  validity  of  these  assumptions  cannot  be  demonstrated  in  any  rigo¬ 
rous  way  on  theoretical  grounds.  What  we  have  done  is  simply  to  make  some 
heuristic  assumptions  concerning  the  possible  behavior  of  solutions  of 
the  Boltzmann  equation  for  wavelengths  long  compared  to  d. 

If  these  qualitative  arguments  are  accepted,  it  follows  that  one  can 

'  t 

estimate  the  loss  of  amplitude  for  a  long  wavelength  wave  upon  reflection 
from  the  transition  layer.  We  are  limited  to  the  case 

yd  <  1  (15) 

i.e.,  where  y  is  the  z  component  of  the  wavenumber:  we  must  have  wave¬ 
lengths  that  are  at  least  2  it  times  the  thickness  of  the  transition  layer. 
Furthermore,  in  order  to  compute  roughly  the  decay  of  the  wave  as  it 
travels  through  this  layer  in  both  directions,  v'  assume  an  attenuation 
like  e  where  L  is  the  effective  path.  For  order  of  magnitude  argu¬ 
ments  we  may  assume,  in  this  part  of  the  frequency- wavenumber  domain, 
moderate  angles  of  incidence  and  x,z  wavenumbers  of  the  same  order.  We 
thus  take  k  =  y  with  a  total  effective  path  length  of  the  order  of  4d, 
i.e. , 

qL  a  4dq  ;  ^  yilyd  (16) 

Thus  the  reflection  coefficient  modulus  would  be: 

K |  e-1l  --  e-4/3  ^  (17) 

3 

Remarkably  enough,  if  we  take  yi  ~  (which  is  consistent  with  equation 
12) ,  we  get: 

1*1  =  e"'Yd  (18) 

This  is  the  result  obtained  by  Vanowitch  (1967)  from  the  continuum 
equations,  for  the  reflection  of  internal  gravity  waves  from  a  region  of 


large  kinematic  viscosity.  Considering  the  relationship  between  kinetic 
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theory  and  continuum  mechanics,  this  is  not  entirely  surprising.  We  have 
simply  confirmed  the  statements  that  most  of  the  attenuation  comes  from  a 
layer  in  which  the  kinematic  viscosity  becomes  prohibit ive,  that  the  thick 
ness  of  this  layer  is  of  the  order  of  a  scale  height  d  and  that  for  wave¬ 
lengths  much  greater  than  d  reflection  takes  place,  lor  very  long  wave¬ 
lengths  ,  (17)  or  (IS)  show  that  this  reflection  becomes  total. 


l-'or  modest  amounts  of  attenuation,  then,  we  write: 

|R|  =  1  -  e 

with 


(19) 


e  =  qL  -  Tiyd  (20) 

If  now  we  consider  a  waveguide  or  thickness  h  in  which  multiple  imperfect 

reflections  of  this  type  will  occur,  then  the  attenuation  in  the  horizon- 

**  6 1* 

tal  direction  goes  like  e  ,  where  (Tolstoy  and  Clay,  1966): 


6  -  e/2h  tane  (21) 
where  it  is  assumed  that  no  losses  other  than  those  due  to  reflection 
at  the  top  surface  take  place. 


It  is  perhaps  more  convenient  to  obtain  an  expression  for  the 
attenuation  in  turms  of  the  wavenumber  and  frequency.  For  this  we  return 
to  (9).  Using,  as  in  (16),  L  =  4d  we  have: 


Thu: 


-  k3 

qL  :  e  -  4n  —  d 
^  u 


.  .  d  k3  1 
6  -  2n  -r-  — 


h  u  tanO 

assuming  modest  0,  i.e.,  tan  6  -  1 ,  we  have: 


(22) 


(23) 


-  d  k3 

5  :  2n  r  — —  (24) 

h  us 

i.e.,  the  attenuation  is  given  by  the  usual  factor,  weighted  by  the  ratio 
of  the  path  length  in  the  highly  viscous  layer  to  the  total  path  length, 
n  is  given  by  (14),  and  it  is  not  easy  to  pick  a  suitable  figure  for  n«. 
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The  graph  given  by  Midgley  and  hicmohn  (1%6)  gives  n  -  i U ,J  m2  sec  1  at 

;  =  200  km.  It  seems,  though,  tiiat  tiie  molecular  viscosity  decreases  with 

1 

height,  somewhat  counteracting  the  1/,  effect.  This  suggests  something 
of  the  order  of  10b  m  2scc  1  at  heights  in  the  400-500  km  range.  On  the 
other  hand,  if  wo  take  an  elementary  kinetic  gas  theory  result  for  a 
small  hard  sphere  model,  we  find  (Tolstoy,  1967)  for  the  mean  free  path 
in  meters : 

l  -  6.2  x  10‘8  p_1  '(25) 

and  (7)  gives 

n  =  j  <v>  x  6.2  x  10  8x  p  1  (26) 

assuming  <v>  -  c  :  8  x  102  m. sec  1  and  p  =10_11kg.  m  3  at  the  base  of 
the  transition  layer,  somewhere  near  the  500  km  level, 

no  -  1.8  x  106  (27) 

Take  then,  in  (24) 

Ho  =  1.4  x  106 

n  -  5.6  x  106  (28) 

Y  -v.15 
h 

and  we  have 

k3 

6-1.7x10  ——  (29) 

u 

Assuming  a  phase  velocity  in  the  neighborhood  of  600  m.sec  1 ,  wavelengths 
of  the  order  of  6  x  102  km.  will  propagate  to  ranges  of  the  order  104  km 
with  an  amplitude  decay  of  .1.  Thus,  with  the  numbers  we  have  used,  only 
wavelengths  in  excess  of  6  x  102  km  will  propaga.te  to  any  distance. 

These  very  simpleminded  calculations  suggest  that  the  waveguide  mode 
picture  for  internal  gravity  waves  could  be  valid  for  long  enough  wave- 
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lengths.  As  is  well  known,  it  also  applies  for  that  part  of  the  short 
wavelength  spectrum  for  which  the  energy  finds  itself  trapped  below  heights 
of  200-300  km  (see  e.g,,  Harkrider,  1964;  Harkridcr  and  Wells,  1968).  It 
seems  therefore  that  a  good  part  of  the  relevant  explosion-generated  inter¬ 
nal  gravity  waves  can  be  so  treated.  However,  it  must  always  be  kept  in 
mind  that,  ultimately,  everything  hinges  on  the  magnitude  of  q:  satellite 

drag  measurements  at  these  heights  suggest  that  this  quantity  will  vary 

» 

a  great  deal,  depending  upon  solar  activity  and  other  factors  (Harris  and 
Priestcr,  1967;  Schilling,  1967).  One  may  thus  expect  the  actual  attenu¬ 
ation  to  vary  substantially  (and  in  both  directions)  from  the  above  esti¬ 
mates,  which  can  only  be  regarded  as  plausible  mean  values. 
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FI  CURL:  CAPTIONS 


Figure  1:  Complete  absorption  of  an  upward  traveling  internal  gravity 
wave  in  an  isothermal  layer  with  kinematic  viscosity  tending 
to  infinity  as  z  -*• 


Figure  2:  Reflection  with  loss  of  amplitude  by  transition  layer  of 
high  viscosity  but  thin  compared  to  the  wavelength. 


